Interplay of internal stresses, electric stresses and surface diffusion in polymer films 
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We investigate two destabilization mechanisms for elastic polymer films and put them into a 
general framework: first, instabilities due to in-plane stress and second due to an externally applied 
electric field normal to the film's free surface. As shown recently, polymer films are often stressed 
due to out-of-equilibrium fabrication processes as e.g. spin coating. Via an Asaro-Tiller-Grinfeld 
mechanism as known from solids, the system can decrease its energy by undulating its surface 
by surface diffusion of polymers and thereby relaxing stresses. On the other hand, application of 
an electric field is widely used experimentally to structure thin films: when the electric Maxwell 
surface stress overcomes surface tension and elastic restoring forces, the system undulates with 
a wavelength determined by the film thickness. We develop a theory taking into account both 
mechanisms simultaneously and discuss their interplay and the effects of the boundary conditions 
both at the substrate and the free surface. 

PACS numbers: 68.60.-p,68.55.-a,83.10.-y 



I. INTRODUCTION 

The stability of polymer thin films is an important re- 
search subject in polymer physics and materials science. 
On the one hand, the aim may be to obtain a stable 
film, as in coatings and lubrification. On the other hand, 
soft films are used for microstructuring, where they are 
destabilized to yield well-designed patterns that are used 
e.g. as a mould for further microfabrication processes. In 
both cases it is crucial to understand the stabilizing and 
destabilizing mechanisms that prevail in polymer films, 
which can be either internal (Van der Waals forces due to 
reduced dimensions, internal stresses, decomposition in 
mixtures) or external (external stresses, external fields). 

In recent studies on spin-coated polymer films it be- 
came apparent that thin films are prone to store residual 
stresses [l|, H| . Such stresses are created due to the fast 
evaporation process of the spin-coating process: as evap- 
oration is fast, the polymer chains do not have the time 
to reach their equilibrium configurations and in the final, 
glassy state the film has frozen-in non-equilibrium con- 
figurations that give rise to stresses. If these stresses are 
not relaxed, e.g. by ageing or tempering the films, they 
influence the film stability as recently shown in dewetting 
experiments and discussed theoretically 2-6]. There it 
has been shown that stresses increase the initial dewet- 
ting velocity and also strongly influence the long time 
dynamics of the dewetting films. In case the film does 
not dewet, the stresses may still lead to destabilization 
0, as they should give rise to an Asaro-Tiller-Grinfeld 
instability @-[Ioj]. This mechanism has been proposed 
for stressed solids in contact with their melt or for solids 
which evolve via surface diffusion. Its origin is the fact 
that the solid can relax stress and lower its energy by 
creating surface undulations. For polymer thin films the 
interplay between residual stresses and other, e.g. exter- 
nally applied, destabilization forces constitutes an inter- 
esting question of importance for all further manipula- 



tions of freshly spin-coated films. 

In this work we reformulate the energy approach usu- 
ally used to describe the Grinfeld instability in a way that 
highlights the connection with other known instabilities 
in thin films. We use the bulk elastostatic equations to- 
gether with a time-dependent kinematic boundary condi- 
tion at the free interface. A direct coupling term between 
the height of the polymer film and the displacement field 
arises, which has not been discussed before as it is less 
relevant in atomic solids. In polymer films, however, this 
coupling should be present and important. Moreover this 
term establishes the connection to other elastic instabili- 
ties, namely to a buckling-likc instability under compres- 
sive stress and, in the case of an externally applied field, 
to the elasto-electric instability investigated by Sharma 
et al. pU4l^ |. Finally the growth rate of the height of 
the polymer film is derived in case of simultaneous ac- 
tion of stress and external field. This result is briefly 
compared to recent experiments concerning the electro- 
hydrodynamic instability of very viscous (high molecular 
weight) spin-coated thin polymer films heated above the 
glass transition . 

The work is organized as follows: First, in section 
HI1 we recall the classical, energy-based formulation of 
the Grinfeld mechanism. In section IIII1 we start from 
a nonlinear elastic theory, derive the bulk elastic equa- 
tions and investigate in section IIVI the stability under 
stretch/compression. In sectionlVlwe show that by allow- 
ing surface diffusion via a kinematic boundary condition 
for the height of the film, the Grinfeld result is regained 
in a well-defined limit. The coupling between height and 
displacement via the kinematic boundary condition can 
influence the classical Grinfeld instability for intermedi- 
ate stresses. In section IVII we add the external electric 
field to our description. We regain the instability dis- 
cussed in Sharma et al. [11] in a certain limit. Moreover 
the full growth rate of the film height is calculated and 
its consequences for experiments are briefly discussed. 
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FIG. 1: Sketch of the geometry. The thickness of the film 
spans from — ho to £(x) in ^-direction. There is uniaxial stress 
in ir-direction, which can be either compressive (as shown) or 
tensile. 



II. GRINFELD INSTABILITY - CLASSICAL 
WAY OF CALCULATION; EFFECTS OF 
BOUNDARY CONDITION AT THE SUBSTATE 

To start with we briefly review the classical treatment 
of the Grinfeld instability of an elastic medium under 
uniaxial stress [13, EH- Usually a semi-infinite solid is 
investigated, but in view of the thin film geometry we 
allow for a finite thickness ho of the film. The known 
results for the semi-infinite case can then be obtained by 
performing the limit ho towards infinity. As the dynam- 
ics of this instability is energy-driven - the system can 
lower its energy by creating surface undulations - all the 
information needed to describe the system is contained 
in the (free) energy of the system, which has an elastic 
part, Ed, and a surface part, E sm { . 

We chose the coordinate system in such a way that the 
free surface is at z = 0, see the sketch of the geometry in 
Fig . [TJ For simplicity we assume a plane strain situation 
[l9( where the uniaxial prestress Co is taken along the 
x-axis. Consequently, we consider an undulation of the 
surface along x given by 



((x) = sAcos(kx) . 



(1) 



We assume either an infinite system or periodic bound- 
ary conditions in x-direction. k is the wave number of 
the perturbation, A its amplitude and e a small book- 
keeping parameter used in the following when dealing 
with expansions. 

The elastic energy of an linearly elastic solid can be 
written via the stress field <j a p as [20| 



^=^/[(l 



'a/3 



va aa <jpp \ dx dz . 



(2) 



with E the elastic or Young's modulus. Summation con- 
vention is implied for indices occurring twice (a, ft — 
1..3). Using a plane strain approximation, one gets 



Ecl = JE J ^ + ^ U%ii ~ VGli<J i^ 



dx dz . 



E 



(3) 



and 



where now i,j = 1..2 (1 «-> x, 2 z) and E = -^—ji 
v = jzt;- Assuming incompressibility, i.e. a Poisson's 



ratio of v = 1/2, one gets E = v = 1. We also will 
use the shear modulus G later on and note the known 



relations, G = 2 (i+ v ) ~ 2(i+v) ^ = As the system 

is invariant in y-direction, E e i has units of energy per 
unit length. 

The second energy in the problem is the surface energy 



E„ 



7 / (yi + C'0*0 2 -l) dx, 



(4) 



where 7 is the surface tension and £surf is measured with 
respect to the state of a flat surface. 

To evaluate the elastic energy, one has to solve the 
elastostatic problem. The prestress is uniaxial along the 
x-axis and given by a Q xx = ctq, a® z = and a xz = cr zx = 0. 
Note that ao < holds for the case of a compressive stress 
and (To > in case of a tensile stress. Undulations of the 
surface will give rise to an additional relaxational stress 
The total stress cry = of^- + cry has to fulfill the 
Cauchy equilibrium equation 



ViOYj = 



and the compatibility equation 

V 2 {a xx + a zz ) = 



(5) 



(6) 



where V = (d x ,d z ). As the prestress cr?- trivially fulfills 
these equations, we introduce the Airy stress function 
\(x , z) for the relaxational stress via the known relations 



d 2 X 
dz 2 



d 2 X 
dx 2 



d 2 x 
dxdz 



(7) 



The equilibrium equation is then automatically fulfilled 
and the compatibility reduces to V 2 V 2 x = 0. This bi- 
harmonic equation has to be solved with the following 
boundary conditions (BC). 

At the free surface z — C,(x), the normal-normal com- 
ponent of stress has to balance the surface tension, while 
the shear stress has to vanish. With n and t denoting 
the unit vectors normal and tangential to the surface, 
respectively, the BC at the free surface read 



C"(x) 

n ^i n i = 7 [ 1 + C / (a .)2] 3 /2 ' = • 



(8) 



or explicitly 



(00 + Vxx)(' 2 - 2&xzC + &zz = 7 

—axzC! 2 + C!{b zz —oo — o xx )+b xz = 0. 



C" 



[i + C' 2 ]V2 



(9) 



Note that all the stresses in Eqs. ([U HJ) have to be eval- 
uated at the interface, i.e. at z = Ci x )- 

At the bottom surface z — —ho, where ho is the film 
thickness, we impose vanishing normal displacement 



u z — at z 



(10) 
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meaning that the film is not allowed to detach from the 
substrate. As the second BC, we study two possibilities, 
depending on the preparation of the system: First, to 
study the case of possible slippage at the lower interface, 
one prescribes 



slip BC 



at z 



-ho 



(11) 



implying vanishing shear stress at the bottom (or equiv- 
alently a vanishing force on the lower surface of the film 
in x-direction, i.e. no traction force). This condition will 
be called 'slip BC in the following. A second relevant 
situation, applying to the case where the polymer film is 
rigidly attached to the lower surface, will be referred to 
as 'fixed BC, 



fixed BC 



at 



-h ■ 



(12) 



We will see that these two different BC, slip vs. fixed, 
have a qualitative influence on the instabilities discussed 
in the following. 

Slip BC at the bottom: The solution of the elasto- 
static problem with the slip BC at the bottom, Eqs. ([TU1 
[TTj) . is the Airy stress function Eq. (|A1|) given in appendix 
IA1 The coefficients occurring therein have to be deter- 
mined by the BC at the free surface: one calculates the 
stresses via Eqs. ([7]), evaluates them at the free surface 
z = C( x ) an d expands in powers of e. From the BC at 
the free surface, Eqs. ((9]), one then determines the coeffi- 
cients in the Airy stress function at order 0(e), yielding 
Eqs. flAjp . 

The problem is now solved at linear order in the un- 
dulation, and we can study the corresponding energy of 
the system. The elastic energy will change due to the 
undulation-induced relaxational stress This change, 
AE C \ = E c \ - E°,, explicitly reads 



AE el = — 
2E 



(TO + &xx) + o\ z + Ad 2 xz 



—2(& xx + a )a zz - (Tq 



dx dz . (13) 



With the Airy stress function determined, the stress field 
can be evaluated. The integrations in Eq. (jTB")) have first 



to be performed over the film thickness, dz. Then 
one usually averages over x, assuming periodic boundary 
conditions: by writing (E) it is understood that one has 

averaged like 57 J dx. Note that due to this averag- 
ing the contribution in 0(e) vanishes. To leading order 
0(e 2 ) one calculates 



(A£el, B > 



A 2 k [a 2 + P7 2 + (a 2 - k 2 ^ 2 ) cosh (2h k)] 



E (2h k + smh(2h k)) 



Note that the averaged quantities, (AE), have units of 
energy per unit area. 

Let us now briefly discuss the obtained result. To 
regain the classical limit of an semi-infinite elastic half 
space one performs the limit ho — > 00. The change in 
total energy, (AE to t) — (AE e \ + E SUI f), then reduces to 



The surface energy is directly calculated from Eq. Q and 
yields in order 0(e 2 ) 



{Esui-i) 



l{^\xf) 



2;„2 



T A Z k 



(15) 



A 2 
' 4 



2 Sk 

E 



-fk' 



(16) 



As becomes apparent, the prestress leads to an decrease 
of the energy. The stress enters quadratically, thus 
both compressive and tensile stress trigger the instabil- 
ity which makes it different from buckling instabilities 
|20| . see also section HVl The surface tension acts against 
the instability and stabilizes high wave numbers, see the 
second term in Eq. (|T^)) . 

As we are predominantly interested in polymer thin 
films, let us perform the opposite limit, hok <C 1. This 
amounts to saying that the modulation wavelength is 
large compared to the film thickness. In this 'thin film' 
limit, the change in total energy reads 



A 2 

(A£ toM ) =- — 



Eh 



1 



\{h k) 2 



jk' 



(17) 



Again, the prestress is destabilizing, independent of its 
sign. However, there is no wave number dependence of 
the destabilizing term to leading order. The same result 
was recently obtained in the framework of a lubrication 
approximation 0. 

Fixed BC at the bottom: In this case the Airy 
stress function given by Eqs. (|A31 IA4I) has to be used. 
For the change in the elastic energy this results in 



(A£kf) 

A 2 k [(a 2 



k 2 j 2 ) sinh (2h k) + 2h k(cr 2 + fc 2 7 2 )] 



E {2h 2 k 2 + cosh{2h k) + 1) 



(18) 



To cross-check, in the classical limit of a semi-infinite 
elastic half space one again recovers Eq. (|Tt)|) . This is 
expected as for a half space the BC at the bottom should 
not be important. In contrast, in the thin film limit 
h k <C 1 one gets 



A 2 fial 
4 \Eh 



(h k) 2 - 7 fc : 



(19) 



Note that the prestress still lowers the energy, but now 
has the same wave number dependence as the contribu- 
tion from surface tension. Thus only above a threshold, 



(14) 



(JO > 




(20) 



the prestress can destabilize the system. 

We have seen that the total energy of the system can be 
lowered by surface undulations in all the cases discussed 
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above. To establish these favorable undulations, it needs 
a mechanism that allows rearrangements to occur. In 
the classical case of a solid in contact with its vapor, 
this is achieved by melting-crystallization processes at 
the surface. This results in a velocity of the boundary 
vmc = [§], where V is a mobility. 

A second possible mechanism - on which we would like 
to focus here in view of polymers - is surface diffusion. If 
atoms or vacancies (in case of a solid) or polymer chains 
(in case of polymer films) feel an inhomogeneous chemical 
potential at the surface, they will diffuse. As a result, 
the boundary will move with a velocity vd = — Mk 2 AE, 
with a mobility coefficient M 8] . Note the second order 
spatial derivative stemming from the diffusion process 
and reflected in the A: 2 -dependence of vd- 

Much more is known about the Asaro-Tiller-Grinfeld- 
instability, for which we refer to the literature. For the 
nonlinear evolution beyond the instability, see [U for 
analytical work and [23] for phase-field modeling. Con- 
cerning experiments, very clean realizations of this in- 
stability have been observed in Helium crystals [24| and 
single crystal polymer films [25| . 

III. STRETCHED ELASTIC SOLID: 
NONLINEAR BULK FORMULATION 

To properly describe finite stresses in a thin polymer 
film, one has to use a nonlinear elasticity formulation. 
Let us assume that the film was originally in a stress- 
free state, described by coordinates X = X;ei. Then we 
stretch (or compress) the film, for simplicity uniaxially 
in the 3>direction by a factor A > 1 (A < 1). This state 
will be described by coordinates x = o^ei and considered 
as the base state. This state will be under uniaxial stress 
a xx , see below. Finally the film is brought in close con- 
tact with the substrate (either still permitting for slip, or 
perfectly fixed to it, see the two BCs discussed in the last 
section) and we let it evolve. This current state will be 
described by coordinates x = a^ei. Note that we discuss 
here only the simple situation where the film is attached 
after the stretch. The situation where the polymers at- 
tach to the substrate while the film is stretched (which 
probably better corresponds to the situation during spin- 
coating) is more involved as the uniaxiality is broken due 
to the presence of the substrate, cf. Ref. [261 ]. 

The total deformation gradient from X to x reads 

„ <9x <9x dx „ „ 

F= = =:F 2 -Fi. (21) 

Here 

Pi =Diag(A,A- 1 / 2 ,A- 1 / 2 ) (22) 

describes the stretching (compression) of the film by a 
factor A > 1 (A < 1). Note that this step must be de- 
scribed in the nonlinear regime, as stresses are finite. The 
second tensor (with I the identity), 

F 2 =I + Vu, (23) 



introduces the usual linear displacement gradient tensor 
Vu = (djUijij in the current state with respect to the 
stretched state. As we are only interested in the stability 
of the base state, here a linearized theory is enough for 
our purposes. As usual we denote with B = F- F T and 
C = F T -F the left and right Cauchy-Green tensors. As 
B is in Eulerian frame we adopt it for the stresses. C 
is in Lagrangian frame and is more convenient for the 
energy definition. Using a Neo-Hookean elastic solid [27j 
, the Cauchy stress tensor is defined as 

a = GB - PI . (24) 

It describes the stress after a deformation in the current 
configuration. P is a Lagrangian multiplier (an effective 
pressure having units of [Pa]) that ensures the incom- 
pressibility condition. In the base state, from Eq. (|2"4")l 
one directly gets cr°^ = except for 

a° xx = G{\ 2 -\- 1 ) . (25) 

This establishes a connection between the stretch factor 
A and the prestress Co- 

Now, let us consider again a plane deformation with 
respect to the prestretched base state. Evaluating the 
Cauchy stress tensor in linear order in the displacement 
gradient, imposing plane strain and using incompressibil- 
ity, one arrives at the bulk equations 

G {\ 2 d 2 x u x + X^dlux) - d x P = 0, 

G {\ 2 d 2 x u z + \~ 1 d 2 u z ) - d z P = 0. (26) 

Note the asymmetry introduced by A =/= 1, i.e. the pre- 
stretch. All quantities can be expressed either in the 
base state x or in the current state x - as deformations 
u are small, they amount to the same expressions. For 
A = 1 one regains the classical elastostatic equation for 
an incompressible solid, GV 2 u + VP = , where u is the 
displacement field. 

The elastic energy density, p c \, for the Neo-Hookean 
clastic solid reads 

AU = | (Tr(C) - 3) . (27) 

Here we did not include the pressure as a Langrangian 
multiplier (giving rise to a term +P(dct(B) — 1)), 
as incompressibility is imposed when solving the bulk 
equations, see the next section. Note that for plane 
strain and small deformations one regains Eq. ([3]) to 
second order in displacement gradients, i.e. p c \ — 

[(1 + v) Oij ~ WiiVjj] ■ 

Now we have established the equations for a nonlinear 
prestretch and a subsequent linear theory. Note, how- 
ever, that the Neo-Hookean model should not be used 
for A- values too far from A — 1. Otherwise effects of 
e.g. the crosslink length must be taken into account and 
one should use more realistic models like the Mooney- 
Rivlin solid [27J . We will now investigate the stability of 
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the prestressed base state with respect to surface undu- 
lations: (i) for the purely elastic case, (ii) in the presence 
of surface diffusion, making a connection with the classi- 
cal Grinfeld instability, (iii) in the presence of an electric 
field normal to the free surface, regaining and general- 
izing results obtained previously [U HH and finally (iv) 
with both surface diffusion and applied electric field. 



IV. STRETCHED ELASTIC SOLID: 
SOLUTIONS FOR SURFACE MODULATIONS 

In this and the following section we solve the elastic 
bulk equations and show how surface diffusion can be 
incorporated within this approach in a generic way to 
regain and generalize the Grinfeld result. We use the 
same boundary conditions as introduced in section [Til 
i.e. Eq. © for the free surface and either the slip BC or 
the fixed BC at the bottom. 

The stability of the base state can be studied by the 
ansatz 



u x (x, z,t) 
u z (x, z, t) 
P(x,z,t) 



u x (z)e lkx+st 
u z (z)e ikx+st 
p{z)e ikx+st 



(28) 



where the amplitudes are small perturbations of order 
0(e) in height perturbations, see Eq. (j3"2")l below. Note 
that we allowed for a temporal dependence which will be 
used only in the following sections. Using incompress- 
ibility, iku x + d z u z = 0, one obtains a single decoupled 
equation for u z given by 



k 4 X 3 u z - fc 2 (A 3 + 1) < + u 



(4) 



0. 



(29) 



With this equation solved, one easily obtains u : 
from incompressibility and the pressure from p(z) = 

(129 



G (y^jjjr- — \ 2 u' z (z)J. The general solution of Eq. 
reads (for A ^ 1 |39f ) 



u z (z) — e ^2 {Ai cosn [h ( z + ho)} 

i=l,2 

+Bi sinh [hi (z + h )] } 



with 



ki = k , and ki = l = fcA 3/2 . 
Imposing the BCs at the substrate yields 



(30) 



(31) 



The normal vector of this surface reads n = 
(—ikh(x, t), 0, 1) at first order. Thus at the free interface, 
cf. Eqs. ©, the BCs read (Jijfij — —■~fk 2 h(x,t)ni. They 
fix the remaining unknown coefficients and one obtains 
A\.s = M,s = and 

2k (G (I 2 - k 2 ) cosh(M) + fc 2 Z7Asinh(M)) 



B hs = -h 



B 



2,s 



G(k 2 + P)g+(k,l) 
fc 3 7 Asinh(/i fc) + G (I 2 - k 2 ) cosh(/i fc) 



lGg+(k,l) 

in case of the slip BC at the bottom; for the fixed BC 
G (I 2 - k 2 ) v(k, I, k) + fc 3 7 Aw(fc, I) 



(33) 



Ai, f = h 



G (fc 4 + 6k 2 l 2 + l 4 - (k 2 + I 2 ) f(k, I)) ' 

G(l 2 -k 2 )w{k,l) +1 \kMk,l,l) 
G (k± + 6k 2 l 2 + I* - (k 2 + I 2 ) f(k,l)) 1 ' 



—B\ r. We introduced the 
I 



and A 2 j = —Aij, B 2 j ■ 
following abbreviations 

5 ± (fc,0 = ±BiDh[(l + k)ho](l-k) 
+ sinh [(/ - fc) ho] (I + k) , 
/(fc,/) = cosh [(/ + fc) ho] (I - fc) 2 

+ cosh[(7-fc)/io] {l + kf , 
v(k,l,m) = (I 2 + fc 2 ) cosh [kh ] - 2m 2 cosh [I ho] , 
w(k, I) = (I 2 + fc 2 ) sinh [kho] - 2kl sinh [Iho] . (35) 

With the general solution obtained, we can now in- 
vestigate whether the base state is stable or unstable. 
According to Eq. (|32p , solutions with non-zero wavenum- 
ber, if they exist, correspond to surface undulations. The 
condition for nontrivial solutions to exist can be written 
as 



h(x, t) = u z (x, z = 0, t) 



(36) 



or eh — u z (0). Namely, for consistency the displacement 
at the surface must equal the height perturbation. An 
alternative formulation would have been to write down 
the system of BCs as a 4 x 4-matrix equation and looking 
for nontrivial solutions via the zeros of the determinant. 
With u z known, Eq. (|36[) can be written as 

hZ(k) = (37) 

with a function of wave number Z(k). If one finds wave 
numbers k* with Z(k*) — 0, periodic solutions exist; oth- 
erwise hZ(k) = implies h — and the film stays flat. 



u z {-ho) 


= 0^A 2 


= -A l 




Explicitly, for the two considered BCs one gets 


&xz(-h ) 


= 0-^A 2 


= 0, 


(slip BC) 


Z s (k) = 4k 3 lG cosh (lh ) sinh (kho) 


u x {-h Q ) 


= 0^B 2 


— jBu 


(fixed BC) . 


-smb. (Iho) (I 2 + k 2 ) 2 G cosh (kh ) 



As before, we parameterize the upper free interface of 
the thin polymer film by a harmonic function with small 
amplitude of order 0(e) 



sinh (Iho) fc 3 7A (I 2 - fc 2 ) sinh (kh ) , (38) 



G 



Z f (k) = 4Gk 2 l(l 2 + k 2 ) + —(r(-l,k)-r(l,k)) 



z = h(x, t) = the 



ikx+st 



(32) 



+ -fc 3 7 A (Z 2 -fc 2 ) g~ (k,l) 



(39) 
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with r (k, I) = (k+l) 2 (k 3 +3kl 2 -k 2 l+l 3 ) cosh [h {k - I)}. 

For both BCs, nontrivial solutions do not exist un- 
der tension, A > 1, as one would expect. Buckling oc- 
curs under compression, but only for non-physical values, 
namely for A < A c ~ 0.03 for a typical surface tension 
of 7 = Q.5hoE. For such high compressions, the Neo- 
Hookean law is no longer a good description. Moreover, 
the assumption that the film stayed flat in the first step 
(from X to x, i.e. before attaching to the substrate) is not 
valid anymore - the film would have buckled long before. 
Indeed the threshold for buckling for two free surfaces 
should be lower than for the BC that the film does not 
detach from the substrate surface, Eq. (ITU1) . Thus we can 
conclude that the film stays flat for all reasonable values 
of A, G and 7. Note, however, that films can be unsta- 
ble if they are swollen in the presence of the substrate, 
cf. Refs. (Mill. 



V. ADDING SURFACE DIFFUSION - THE 
GRINFELD INSTABILITY AGAIN AND 
CORRECTIONS 

In the last section we investigated the stability of the 
base state with respect to in-plane stresses and found 
that the purely elastic system is stable. Here we add 
the effects of diffusion of polymer chains close to the film 
surface due to stress relaxation-induced changes in the 
chemical potential. As a consequence the system can 
produce undulations by diffusive transport of material, 
in addition to possible elastic displacements. We show 
that one regains the Grinfeld instability in a well-defined 
limit. The overall result is more general as it comprises 
corrections to the Grinfeld mechanism, see below. 

If we allow for surface diffusion, Eq. (|36|) has to be 
modified in order to allow for this dynamics. For the 
height modulation h = h(x,t) one can write 



d t h = d t u z \ z=0 - (d t u xlz=0 )(d x h) 
+Md 2 x (S^ z=0 ) . 



(40) 



The first two terms on the r.h.s. stem from the stan- 
dard kinematic BC at a free surface, usually written as 
dth = v z — v x d x h with h the height of the surface and 
(v X iV z ) = d t {u x ,u z ) the fluid velocity [3(|. The second 
term is purely nonlinear and can be neglected in the fol- 
lowing linear analysis. 

The last term on the r.h.s. represents the surface diffu- 
sion (note that in three dimensions d 2 has to be replaced 
by the surface Laplacian (ll[). ^ wm smoothen gradients 
in the chemical potential, which is given by 



Sfi = SE e 



(41) 



k is the mean curvature of the surface, given at 0(e) by 
k = d 2 h. SE e \ is the change in elastic energy density due 
to the surface undulation, compared to the flat surface. 
The coefficient M is a mobility [32|, |33| and explicitly 

reads M = D ^ s ^ , where fcgT is the thermal energy, 



V is a microscopic volume (of the polymer chain in our 
case), D is the surface diffusion coefficient and n s is the 
surface density of diffusing objects. Note that in the 
view of recent experiments on spin-cast polymer melts, 
we here allow for a finite chain mobility (at least close to 
the free surface), although we assumed a purely elastic 
behavior of the film. A generalization of our approach to 
the more adequate viscoelastic case will be the subject of 
a forthcoming study. 

Eq. (|4"01) for the dynamics of the surface undulation 
is further motivated in appendix [Bj The terms arising 
naturally from the kinematic BC are commonly not in- 
cluded in the treatment of the Grinfeld instability, as in 
the usual context one concentrates on the diffusive trans- 
port of atoms or vacancies. Taking the coupling to the 
displacement into account - if extended objects like poly- 
mers are diffusing - corrections to the 'classical' Grinfeld 
behavior arise: the time derivative in dtu z \ z=0 leads to 
a renormalization of the growth rate s(k) of the height 
perturbations h(x,t) = ehe lkx+st . In view of this, in the 
following we will sometimes compare the 'classical' Grin- 
feld and the 'kinematic' case. 

In the previous section we have already calculated the 
general solution for the displacements. Thus the stress 
tensor is also known and using Eq. (|27l) one gets the 
changes in the elastic energy 5E e \ = p c \ — with respect 
to the base state 



SE„ 



SE, 



eG (I 2 - k 2 ) cos(kx) 
■ s = PA 

(Si, s fc cosh [k {h + zj\ + B 2 , s l cosh [I (ho + «)]), (42) 
eG (I 2 - k 2 ) cos(kx) 



elj 



k 2 \ 



■ (Aij (fcsinh [k(ho + z)} - Ismh [l(h + z)]) 
+B 1 . f k (cosh [k(h + z)] - cosh [I (ho + z)]) ) ■ (43) 

For the surface energy, as before Eq. (g]) yields E sur f = 
ecos (kx) e st ^/k 2 h. Now we can proceed in two ways: 

Classical calculation, nonlinear case: First we 
can use the classical Grinfeld argument, i.e. we inte- 
grate from — ho to h(x) over the film thickness and av- 
erage over the assumed periodic x-direction to obtain 
A-Ed = (J 5E C \ dz) x . Upon averaging the linear order 
in e vanishes. At 0(e 2 ), one gets to leading order in k 



SE e 



h 2 e 2 (A 3 — l) (3A 3 + l) G 



SE, 



16h X(X 3 + if 
h 2 k 2 e 2 h (A 3 -1) 2 G 
4A 



0(k ) , (44) 



+ 0(k 4 ) . (45) 



Let us compare to the result obtained in section HH In 
the limit A = 1 ± 6 with S <g^l and using 8 = <r /(3G) as 
implied by Eq. (|25[) in this limit, one gets including the 
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surface energy 



5E< 



tot,s 



"6 



SR 



tot,/ 



4 \Rh 

h 2 ( w 

' 4 



i + (hoky 



(46) 
(47) 



Note that in leading order this is exactly Eqs. ([171 IT9]) , 
The correction (/tyfe) C 1 in Eq. (|4"6")l has a slightly 
different prefactor as in Eq. (fT7| , which is due to the fact 
that the fully linear calculation from section [TT| is only 
correct for infinitesimal stresses. 

Consistent calculation at order 0(e): The use of 
an averaging in the Grinfeld calculation seems not nec- 
essary to us. We will thus determine the growth rate of 
surface undulations by using 



d t h = dtu z \ z=0 + Md 2 (Sfi\ z=0 ) 



(48) 



The l.h.s. and the first term on the r.h.s. are of first 
order in e. Thus it is sufficient to determine the change 
of the chemical potential at this order, i.e. evaluating Sfi 
at the surface. The full growth rates obtained by this 
equation are given by Eqs. (|A5|) - (|A8|) in appendix lAl for 
the slip and the fixed BC, respectively. In the thin film 
limit hok <C 1, one obtains 



s s (k) ~ Mk 
+Mk 4 



2 G(X 3 l) 2 



2h X(l + A 3 ) 
Gho(X 3 - If 



6 A 



-7 1 



A 3 -l 
2 (1 + A 3 ) 



(49) 



for the slip BC at the bottom and neglecting the 6\u z | z= o 
term at the free surface. Including the kinematic term 
yields 



Ss^kin (^) 

+Mk 4 



~ Mk^ G(A3 ' 1)2 
h X(3 + A 3 ) 

G7i (4 + 3A 3 + A 6 )(A 3 - l) 2 

3A(3 + A 3 ) 2 



- 47 



(1 + A 3 ) 



3\21 



For the fixed BC we get in both cases 



Mk 4 



^(A 3 -l) 2 - 7 
A 



(3 + A 3 ) 2 

(50) 



(51) 



Let us first discuss the limit of small stresses asgain, 

A = 1 ± S with S = cr /(3G) < 1. Both Eqs. (gH [50]) 

2 

yield at leading order in the stress s(k) ~ Mk 2 -j^. Ex- 
cept for a factor of 4, at leading order this is exactly 
Eq. HT7J). The same is true for the fixed BC and Eq. ([T9"]l. 
Hence in the low-stress and low-wave number limit, our 
results obtained for the dynamic BC at the free surface 
are identical to those obtained by the energy-based cal- 
culation in section [TT] in the following sense: s(k) are 
growth rates as calculated from a dynamical equation for 
the surface undulation. When comparing to the Grinfeld 
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FIG. 2: Growth rates as a function of reduced wave num- 
ber for different stretch factors A, corresponding to different 
prestress. Panel a) displays the case of the slip BC at the 
bottom. Without stress, A = 1, the system is stable. Any 
compression (A < 1) or stretch (A > 1) lead to an instability, 
but with differing rates. A = 1.1 and A = 0.9 correspond both 
to a prestress |cro|/S ~ 0.075. Solid lines have been obtained 
with the kinematic BC at the free surface, pointed lines just 
with the surface diffusion. Panel b) shows the case of the 
fixed BC at the bottom. An instability only occurs beyond 
critical A-values. Parameters: 7 = 0.5hoE. 



calculation, there too one has to impose a diffusion dy- 
namics driven by the decrease in energy. One can write 
d t A = -Mk 2 J^(AE tot ), with E tot oc A 2 . The varia- 
tional derivative with respect to A yields a factor of 2. 
Taking into account that in the energy approach one has 
averaged over (cos 2 (fca;)) yields another factor of 2, which 
explains the differing prefactors. However, one should 
note that using the spatial averaging process implies a 
calculation order 0(e 2 ), while our method is 0(e). 

Let us now discuss the effect of finite stretches and 
of the kinematic contribution. At leading order in the 
stress, both Eqs. (|49l [50l) reduce to the Grinfeld re- 
sult. However, in next order in the stress Eq. (|49[) 

yields — \Mk 2 , while the kinematic version yields 



- 77 Klk 2 f^Jl — . 

3 E 2 h 



E 2 h ' 

First, this shows that the symmetry with 
respect to the sign of the stress, i.e. whether it is due 
to stretch (er > 0) or compression (oq < 0), is broken 
by the elastic nonlinearity. Second, the sign of the cor- 
rection is sensitive to whether the kinematic BC at the 
free surface is important (e.g. for diffusion of extended 
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FIG. 3: The fastest growing wavenumber kmax as given by 
the maximum of the growth rate s Si fe in shown in Fig. [2k.) as a 
function of rescaled prestress \oq\/E. Solid lines are obtained 
for 7 = 0.7/to-E, dashed lines for 7 = O.hhoE. 




FIG. 4: Sketch of the capacitor geometry. There is uniaxial 
stress in x-direction (which could also be tensile). An ex- 
ternal electric field in z-direction (for the unperturbed film) 
is imposed by externally applying a voltage V between the 
electrodes of distance d. 



residual stress and an electric field, both acting as desta- 
bilizing factors for elastic films. 



objects like polymers in a network) or not. 

To compare to a real system, we use the following pa- 
rameter values as suggested by Ref. [15]: ho — 140nm 
for the thickness of the film and E = 5 • 10 5 Pa for the 
modulus. For the surface tension we use the value for 
polystyrene, -fps — 30 ■ 10 -3 Nm -1 . Fig. [5] displays the 
full growth rates, Eqs. (| A5|) - (| A8|) in appendix [M as a 
function of reduced wave number khg. Note that we 
renormalized s(k) — (h^/(ME))s(k). Fig. UK) displays 
the case of the slip BC at the bottom, with (solid curves) 
and without (dotted curves) accounting for the kinematic 
BC at the free surface. Finite stresses lead to a Grin- 
feld instability. Growth rates differ whether compression 
(A < 0) or extension (A > 0) is considered. In case of 
the fixed BC at the bottom, see Fig. [5Jd), there exists a 
threshold stress beyond which the system becomes un- 
stable. For the chosen surface tension, 7 = Q.hh^E in re- 
duced units, the system destabilizes for A > A^ c ~ 1.387 
and A < A2, c — 0.426. Note that a (symmetric) thresh- 
old stress also occurred in the linear model, cf. Eq. (|20[). 
Fig. [3] displays the dependence of the fastest growing 
wavenumber on the prestress \ao\/E as obtained from A 
by Eq. (|25[) . One clearly sees the asymmetry with respect 
to compression/stretch for finite stresses. 

To summarize, in the last two sections we proposed a 
general framework that includes the Grinfeld instability 
as well as possible buckling. The possibility of buckling 
is due to the coupling of surface undulations and the dis- 
placement field via a kinemtic BC at the free surface. 
One gets corrections to the Grinfeld instability, as con- 
tained in the full growth rates given in appendix[X] How- 
ever, in the small wave-number limit and for thin films, 
the leading order terms are identical with the classical 
result. For finite stresses the zb-symmetry with respect 
to stresses predicetd by the linear Grinfeld-theory is no 
longer valid. In the next section we use the developed 
framework to study the simultaneous action of in-plane 



VI. ADDITION OF EXTERNAL ELECTRIC 
FIELD 

Recently the instability of polymeric liquids [HI, [H[ 
and elastomers [lij in an external electric field acting 
normal to the film surface has been investigated experi- 
mentally. In Ref. 15j, it has been found that the insta- 



bility is faster for freshly spin-casted films than for aged 
films. This suggests that stresses in the fresh films due to 
the noncquilibrium production process may be involved 
in the destabilization. In view of this we generalize the 
developed approach to the case where an external electric 
field is acting normal to the surface, in addition to the 
stress in x-direction. The electrostatic part will be closely 
related to previous studies of elastic instabilities [ll|-[lj] 
due to forces normal to the surface (Van der Waals or 
electric field). Related studies have been undertaken in 
Refs. [34], HI!- However, there the thin film was regarded 
as conductive, the external stress was imposed externally 
(implying that the base state with applied field was fixed 
at a xx = (To rather than a® x = <jq + F as in our case with 
F the additional contribution from the electric field, see 
below) and the kinematic BC (i.e. the coupling of film 
height and displacement field) at the free surface was not 
taken into account. 

Let us assume that the polymer film is brought into 
a parallel plate capacitor, see the sketch in Fig. [4] A 
voltage difference V is applied over the distance of the 
two plates d (the lower plate is at z = — ho, the upper one 
at z — d — ho). The gap may be filled with any dielectric. 
In view of the experiments in Ref. [llj], we take e\ ~ 2.5 
(polystyrene) as the dielectric constant of the polymer 
film and 62 = 1, i.e. the gap is filled with air. 

The electric field will introduce a stress at the polymer- 
air interface. Let us introduce the Maxwell stress tensor 



(52) 
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where the index k = 1 denotes the polymer film and 
k = 2 the gap. The BC at the free surface, cf. Eq. ©, 
now reads 

n • (a« - tr&> + T« - T< 2 )) • n = jd xx h , (53) 

where we wrote only the linear order expression for the 
surface tension. We can put — (or to a constant 
pressure value that is not important), = a and de- 
fine an electrostatic 'pressure' (strictly speaking a normal 
stress) by 



PE (h) = h- (r«-T( 2 )) -n. 



(54) 



E 



, Fixed BCs 




Note that this electrostatic stress depends on the film 
thickness, see below. The BC finally reads n • a • n = 
<yd 2 h-p E (h). 

We now have to evaluate the additional contribution 
from the electric field. We can again solve the problem 



by a perturbative method by writing = E^ + E^, 

(i) 

where E^ is the first order correction due to undulations. 
To lowest order, one has to satisfy that the normal dielec- 

i!q z . Second, 
= and thus 



(0 



trie displacement is continuous, tiE^l = E { 



we have V = hoE^l + (d - h Q )E^' z , E { 



(2) 



0,x 



one gets E$ = e 2 V/(h + (d - /*>) d), E ( 0z - -^ z 
In the next order, we have to solve Maxwell's equations 



h(2) _ 



(i) 



d z E\ 



(i) 



d z E\ 



■B g E®=0 : 



(55) 



with the BCs 

'1,: 

n- (e 2 E( 2 



E\ l(z = -ho) 



eiE 



= 
= 



E^l(z = d-h )=0, 
t-(E< 2 >-E w ) = 0.(56) 



These BC state that the field has to be perpendicular 
to the conductive electrodes and that at the film surface 
one has continuity in the normal displacement and the 
tangential field. Assuming e[ z cx cos(fex), the system 
is readily solved yielding the field components given by 
Eqs. (|A10I) in appendix [3] in agreement with Ref. 36!]. 

Evaluating the normal-normal component of the 
Maxwell stress, for the electrostatic pressure as defined 
in Eq. ([54| above we get to leading order 



■p E (C) - F + Y(k)(, 



(57) 



FIG. 5: Stability diagram for the electric field-induced insta- 
bility. The curves display the electric contribution to stress in 
reduced units, Yoho/E, as a function of reduced wave number 
hok for the slip and the fixed BC at the bottom, as indicated. 
Solid lines correspond to the stress-free case, A = 1.0. Dashed 
lines are for finite stretch, A = 1.1, and dotted lines for finite 
compression A = 0.9. Parameters: 7 = 0.5hoE. 



an isotropic pressure [lllQjt, given by F = — pe(C — 0). 
Concerning the displacements relative to the base state, 
the procedure is completely analogous to the one in the 
previous sections. Only the BC at the free surface, and 
the chemical potential have to be changed accordingly to 
include the electric stresses. In the chemical potential, 
Eq. (|4"T1) . we have to add the contribution due to the 
electric stress by writing 



5(i = 5E e i - 7« +Pe(C) 

= SEci+j^C-F -Y(k)(. 



(60) 



As F is a constant, its contribution to surface diffusion 
vanishes. 

The general solution for the displacement field, 
Eq. (|3TJ)) , with the BCs at the substrate already imposed, 
is still valid. One only has to determine the coefficients 
fulfilling the new BC at the free interface. These coef- 
ficients, Bi s , Bf s and Af *, Bf* can be obtained from 
the respective solutions without field by the simple sub- 
stitution 



Y(k) 
k 2 



(61) 



where 



F 

Y{k) 



z (e 2 /io + {d- h )ei) 
-2fcp £ (0)(ei -6 2 ) 
[ei tanh ((d — ho) k) + £2 tanh (hok)] 



Note that both F and Y are strictly positive, F, Y > 0. 
As one has n • (T^ - T^) • t = 0, the tangential BC at 
the free surface is unchanged by the electric field. In the 
base state, the contribution of the electric field will be 



This rescaling of 7 permits to obtain also the functions 
Zf(k), Z^(k) that determine the stability of the flat base 
(58) state in the presence of a field, as well as the growth rates 
sf and . The obtained expressions are very general. 
Although unsightly they contain the physics of buckling, 
the elasto-electric instability, the Grinfcld-instability and 
surface diffusion in an applied electric field. 

Let us first discuss the case without surface diffusion. 
One expects to get an instability for Y > Y c as described 
by Sharma et al. [11]. Note, however, that in case of 
an applied electric field Y(k) is fc-dependent, while in 



(59) 
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Ref. [TT| a Van der Waals-interaction with a contactor 
was studied, where Y is a constant. To compare we write 
Y (k) as a function of Yq where Yq = lim Y (k), 



In the limit of small stresses, A = 1 ± S with 6 = 
<7o/(3G) <C 1, to lowest order k 2 and up to third order in 
stress one gets 



Y(k) 



kY (dei + h fa - ex)) 



ei tanh(fc(c2 — ho)) + £2 tanh(Ziofc) 
One gets the following conditions for instabilities 



(62) 



Y 



Yn = 



—Z s (k) fa tanh [k (d — ho)] + e 2 tanh [fe/io]) 
Xk 2 (I 2 - k 2 ) sinh [kh ] sinh [lh ] fad + fa - e 2 ) ho) ' 

(63) 

(fc) (ei tanh [k (d — ho)] + e 2 tanh [kho]) 



Xk 2 (I 2 - k 2 ) g- (fc, fad + fa - e a ) ho) 



for the slip and the fixed BC, respectively. 

Fig. [5] represents stability diagrams for the elasto- 
electric instability, for both BCs as given by Eqs. (|B"3")) . 

. For high enough Y , there exist solutions with finite 
k. However, as Y(k) depends on k and one has a compli- 
cated dependence on both A and 7, we could not obtain 
simple formulas for the threshold. For Y (k) = const (as 
for Van der Waals-interactions) we find the same result 
as given by Ref. [37| (slip BC) and as given in Ref. QUI] 
(fixed BC). We can observe the following general trends 
due to finite stretches: considering the small wave num- 
ber branch, to get an undulation with the same small 
wave number, Y (X = 0.9) > Y (X = 1.0) > Y (X = 1.1). 
Thus compression acts stabilizing and tension destabi- 
lizing on small wave numbers. On the other hand, for 
the large wave number branch one has Yq(X = 0.9) < 
Yo(X = 1.0) < Yo(X = 1.1), thus tension is stabilizing 
and compression destabilizing. 

Let us now look at the case with surface diffusion. In 
the thin film limit hok <C 1 one gets for the slip BC 



S f(k) ~ Ss (k) + Affc 2 r |^±il 

Y h X 3 /3 



-Mk 4 



2(A 3 + l) + (3A 3 + l)y 2 : 
2 



^ G G(A3-1) +W(1 + 3A2) 

s sMn {k) _ Mk — G(3 + A 3 ) _ /lo y oA 



(65) 
(66) 



excluding and including the effects of the kinematic BC, 
respectively. In the latter case we show only the leading 
order contribution in k 2 . For the fixed BC one gets 



sf(k) ~ s f (k) +Mk 2 Y 
hi (A 3 -l) 



+Mk 4 Y 



+ Y 2 \, 



(67) 



sf Mn (k) ~ s fMn (k)+Mk 2 Y 



MX 



+Mk A ( Yohl (A 3 - 1) + Y 2 ^ + Y 2 ) , (68) 



where we introduced Y 2 



\^Y{k\ k= o. 



sf(k) 



s s,kin(k) 



Mk' 



v 1 Y ° 

^o + co-gr 



, 3E - 2Yr 



3E 2 h Q 



2cr; 



. 9E + 10F 

27Ehn 



(69) 



In 



6E - 5Yc 







3(^-F )^ 



(64) 



,9£ 3 - 12S 2 y - 12EY + 16Y 

9Eho{E-Y ) 3 
,9E 4 - 76E 3 Y + 116E 2 Y 2 - MEY 3 + 16Y 



27E 2 h (E - Yo) 



(70) 



where we introduced Yn — hnYo, and 



\k) 



(k) 



Mk 2 Y . 



(71) 



For the slip BC, Eqs. (1551 16"6"|) display a coupling be- 
tween prestress and the applied electric stress at the or- 
der O (k 2 ). Hence the application of the field breaks the 
icg-symmetry already in lowest order in stress. In the 
small stress limit this coupling is linear like Yo^o- As one 
usually has h Y E (otherwise the elasto-electric in- 
stability takes over), for small wavelengths one might be 
driven to the conclusion that the coupling between the 
electric field and the stress like YoCo implies that com- 
pression acts stabilizing while stretch acts destabilizing. 
However, the destabilizing contribution from the electric 
field, Mk 2 Yo, is usually dominating and thus s E > s, 
compare e.g. Figs. [2£i) and|6ji). Thus the influence of 
the coupling is observed rather beyond the maximum 
of the growth rate, cf. Fig. [S^), as a gap between the 
curves with compression and stretch. However, in case of 
Co <§C hnY <C E, the two destabilizing forces do not add 
and the growth rate under compression is indead slightly 
smaller than the growth rate of an unstressed film. As 
in the case without field one observes finite stress effects 
yielding positive contributions like +ctq with the kinetic 
BC and negative ones like — Oq with the non-kinetic ver- 
sion. Positive and negative prestress have thus opposite 
effects and these effects depend on the kinematic BC. 

Fig. EJa,) displays the general growth rates for the slip 
BC. For A < 1 (do < 0), the kinematic version yields 
a smaller growth rate than the non-kinematic version. 
However, for A > 1 (00 > 0) the opposite is true. Fig. [6)d) 
shows that for the fixed BC finite stress only leads to 
higher order corrections, since the leading order destabi- 
lization is +k 2 Yo while the stress contributions are pro- 
portional to fc 4 . Visible differences between Sf^in (k) and 
Sf (k) appear only for rather high stretch factors, namely 
A > 1.5 or A < 0.7 for the chosen surface tension. For 
Fig. [5]we used again parameters as suggested by Ref. [lj| , 
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a) 




FIG. 6: Growth rates as a function of reduced wave num- 
ber for different stretch factors A, corresponding to differ- 
ent prestress, and with finite voltage applied normal to the 
film. Panel a) displays the case of the slip BC at the bottom. 
Without stress, A = 1, the system is slightly unstable due 
to field-induced diffusion. Compression (A < 1) or stretch 
(A > 1) lead to a more pronounced instability, but with dif- 
fering rates. Solid lines have been obtained with the kine- 
matic BC at the free surface, pointed only with surface diffu- 
sion. Panel b) shows the case of the fixed BC at the bottom, 
where finite stresses only lead to small corrections. Parame- 
ters: 7 = 0.5h E, ei = 2.5, e 2 = 1, d/ho = 5, h Y = 0.001S. 
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TABLE I: Summary of the leading order destabilization terms 
in the growth rate of surface undulations, s(k), for different 
BCs at the bottom and with or without electric field. Both 
stress and electric field are assumed small, <tq, hoY <C E. 



coating [38], etc. As a next step we plan to generalize 
the approach proposed here to the viscoelastic case, to 
come closer to these experiments. 

Another interesting point is that the compressive- 
tensile symmetry holding for the stress in case of the 
Grinfcld instability is broken in several ways: i) by finite 
stresses, ii) due to the kinematic BC, i.e. the coupling 
of film height and displacement field, and iii) due to the 
presence of the external electric field. While the effect 
is of order Oq in the absence of an electric field, it is of 
order YqUq in the presence of field. Thus especially in an 
external electric field, surface undulations have notice- 
ably different growth rates and this may be used experi- 
mentally to determine whether stresses in thin films are 
compressive or tensile. 



VII. CONCLUSIONS AND PERSPECTIVE 



namely a electrode distance of d = bho, dielectric con- 
stants £2 = 2.5 and e% = 1 for the PS film and the air 
gap, respectively, and a voltage of V = 16V. In reduced 
units this leads to 7 ~ 0.5h E, Y h ~ 0.0013i£. 

Let us briefly discuss the relation of this work to 
the experiments of Ref. [HI]. There it has been found 
that freshly produced films, that are supposedly stressed 
due to the nonequilibrium preparation process of spin- 
coating, have faster growth rates than aged films - which 
had time to relax residual stresses. This is in accordance 
with our findings that the two destabilization mecha- 
nisms, the Grinfeld mechanism and the electric force act- 
ing on the free surface of the film, in general join forces. 
However, in Ref. [15| it has been found that the wave 
number of the instability is smaller for fresh films than 
for aged films. This is in contrast to our calculations, as 
in the general case the unstable wave numbers increase 
with stress, see also Fig. [51 There are several possible 
reasons for this discrepancy, like viscoelastic effects in 
the film, inhomogeneities, crust formation due to spin- 



We have studied the instability of a polymer film un- 
der the simultaneous action of internal stress and an ex- 
ternally applied electric field. For this purpose we for- 
mulated a general framework that has a very rich phe- 
nomenology: in absence of surface diffusion the system 
is stable against buckling but displays an electrically in- 
duced instability towards periodic undulations. In case 
that the polymer chains are able to diffuse close to the 
surface due to gradients in the chemical potential, the 
Grinfcld mechanism becomes active, as well as a destabi- 
lizing contribution induced by the external electric field. 
The growth rates of surface undulations are sensitive to 
the boundary conditions at the bottom and have a rich 
phenomenology, see Table |U 

Our approach also highlights the importance of the 
coupling between the height of the film's surface and 
the displacement field inside the film, which naturally 
arises from the kinematic boundary condition at the film 
surface. This coupling has been neglected in previous 
studies. Its consequences can be seen as finite stress 
corrections to the Grinfeld instability, and analogously 
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for the electric instability. Moreover, this coupling es- 
tablishes the connection between the above mentioned 
elasto-electric instabilities and the Grinfeld-like diffu- 
sive instabilities, as becomes apparent from the general 
growth rates of height fluctuations. These growth rates 
have been calculated as a function of internal stress, elec- 
tric field, mobility of the chains and surface tension. It is 
shown that both destabilizing factors, internal stress and 
electric field, generally add. 

The relevance for recent experiments on spin-cast thin 
polymer films has been only briefly discussed. A general- 
ization to the viscoelastic case, and possibly also includ- 
ing more structural details of spin-cast film, is needed to 
account for these experiments. In turn, as the experi- 
ments can measure separately the most unstable wave- 
length and the growth rate, they could give direct access 
to the internal stress and to the mobility of polymers in 
thin films, which both are of technological importance. In 
particular, as the electric field makes the breakage of the 
compressive-tensile symmetry of the Grinfeld instability 
induced by the coupling to the displacement field notice- 
able, careful measurements of the growth rates could be 
used to determine the nature of the stresses, i.e. whether 
they are compressive or tensile. 

The authors would like to thank Ken Sekimoto for 
stimulating discussions. 



Appendix A: Results for stress functions, 
displacement coefficients, electric field 

For the Airy stress functions in the Grinfeld calculation 
(cf. section |n]) one gets 

Xs(x, z) — e cos(kx) [ a s cosh(fc(z + ho)) 

+b s (z + h ) sinh(fc(z + ho))] (Al) 



in case of the slip BC, with 



= 2A 



-2A 



/io/c7cosh (h k) + (7 + uah ) sinh (hok) 
2hok + sinh (2hok) 
kj sinh (hok) + ctq cosh (hok) 



2h Q k + sinh (2h Q k) 



(A2) 



In case of the fixed BC, one calculates 



Xf(x, z) — e cos(kx) [a,f cosh(fc(z + ho)) 

+bf(z + h ) cosh(fc(z + h )) 
+a } k{z + ho)e- Hho+z) ] (A3) 



ith 



a f = 2 A 



-2A 



{ho<Jo + 7) cosh(/iofc) + hokj sinh(hok) 
2/i 2 fc 2 + cosh(2/i fc) -I- 1 ' 
(<to + kj) cosh(/i fc) + h k(a - kj)e~ h " k 



2h 2 k 2 + cosh(2h k) + 1 



(A4) 



The full growth rates calculated in section [V] read 



s s (k) 

Ss.kin (&) 

Sf (k) 

(k) 



M 



2G (I 



' cosh [h k] cosh [hoi] - k 3 jX [b (k, l) + b (k, -I)] 



AUG 



2G (I 



2X (k 2 + l 2 )g+ (k,l) 

2 u2\ 3 



cosh [h k] cosh [hoi] - jk^X [b (fc, I) + b (fc, —I)] 



-M 



X [2 7 fc 3 A (I 2 - k 2 ) sinh [hok] sinh [hoi] + IG (3k 2 + I 2 ) J x (k, I) + kG (k 2 + Wj g~ (k, I)] ' 
(8kH 2 ~fX (k 2 + I 2 ) - k 3 ~/X [r (k, l) + r (fc, -I)] + G (I - k) 3 (I + kf g (k, /)) 
2A (fc 4 + 6k 2 l 2 + Z 4 - (k 2 + I 2 ) f (k, I)) ' 
(8fc 4 Z 2 7 A (fc 2 + I 2 ) - fc 3 7 A [r (fc, l) + r (fc, -/)] +G(l— fc) 3 (/ + fc) 3 gi (fc, I) 



-MIG± 



X [(I 2 - fc 2 ) A 7 fc 3 £r (fc, l) + G(r (-1, k)-r (I, fc)) + 8GZfc 2 (I 2 + fc 2 )] 



(A5) 
(A6) 

(A7) 

(A8) 



where we have introduced linear order. From the Maxwell equations with suitable 

b(k, I) = (fc 4 + 2fc 2 Z 2 + 4/d 3 + I 4 ) sinh [(/ - fc) ho] ■ (A9) 
In section IVII the electric field has to be calculated to 
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BC, as given in the main text, one gets 



(i) _ 



E 



(i) _ 

1,3 _ 



E 



(2) 

1.2 



E (2) 



~ cosh(fc(z + /i )) 

ehcos(kx)e 2 E , 

cosh (/lQ«) 

». • ci \ pSinh(fc(z + /i )) 
— eh sm{kx)e 2 E ■ 



cosh (/lofe) 

~ cosh (fc (z + /i -d)) 
efrcos(M^ cosh(fc(d , M) , 

~ sinh (fc (z + ho — d)) 



-eh sin(kx)eiE 



cosh (k(d — /iq)) 



,(A10) 



x x+dx 
solid at time t=0 




u x \x,t) 



H.,.(.r -Id.r.f) 



x x+dx 
solid at time t 




solid at time t+dt 



with the abbreviation 



FIG. 7: Sketch for motivation for kinematic BC. 



E 



V 



(e 2 h + (d- h ) ei) 

k(ei - e 2 ) 



e\ tanh(/c(d — ho)) + e 2 tanh(/iofc) 



(All) 



at first order 0(e) (which allowed us to replace the up- 
per integral bound h (L(z, t), t) by h (x, t)). The volume 
entering by diffusion along the free surface on L[z = 0, t) 



Appendix B: Alternative motivation for the 
kinematic BC 



Here we want to give a more explicit motivation for 
the kinematic BC with surface diffusion. Let us consider 
a small part of an elastic material bounded by x and 
x + da; on the x-axis, h (x, 0) = and — ho on the z- 
axis, see Fig. [7^,). The extension in y-direction, L y , is 
assumed not to change in time, to stay within the plane 
strain situation. At time t, the respective bounds are 
between x + u x [x, z, t) and x + dx + u x (x + dx, z, t), as 
well as h(x + u x (x, 0, t) , 0) and — ho, see Fig. Let 
us use this state, excluding the upper free surface, as a 
control state and calculate its evolution in time. 

The volume inside the considered piece of material is 

V(t) = h{x + u x {x,0,t),t)L y 

[x + dx + u x (x + dx, z, t) — x — u x (x, z, t)) . 
V{t + dt) = h(x + u x (x, 0, t + dt) ,t + dt) L y 

(x + dx + u x (x + dx, z, t) — x — u x {x, z, t)) 

at times t and t + dt, respectively. The volume change 
SV = V(t + dt) - V(t) thus reads 

SV = d t (h + u x (x, 0, t) d x h) dtL y dx (1 + d x u x ) 

w dthdtLydx, (Bl) 

at first order 0(e). As h(x,t) can be any function, this 
is true also if there are additional surface processes. 

On the other hand, we can express the change in vol- 
ume by changes due to displacements and surface diffu- 
sion on both boundaries. During an infinitesimal time dt, 
the volume change due displacements on the left bound- 
ary L{z, t) — x + u x (x, z, t) is 

ph(x,t) 

— / [u x (x, z,t + dt) — u x (x, z,t)]dzL y (B2) 



(x + u x (x, t) , t) dtLy ~ J S uxf,x {x, t) dtLy. (B3) 

During the same time dt, at the right boundary R(z, t) = 
x + dx + u x (x + dx, z, t) we get 



[u x (x + dx, z,t + dt) — u x (x + dx, z, t)] dzL y , (B4) 



where again we simplified the integral bound, 
h(R(z,t),t) to h(x + dx,t). Surface diffusion at 
the right boundary contributes 

"^surf,a; (x + dx + U x [x + dx, t) , t) dtL y 

— J S uTf,x (x + dx, t) dtLy . (B5) 
From SV — (change at R) — (change at L) we obtain 
8V w —d x J snx f tX dtdxL y 

/ Mx,t) \ 

-dx L y dtd t \ I d x u x (x, z, t) dz . (B6) 



Using the incompresibility condition, integrated over the 
film thickness, one gets 



h(x,t) 



d x u x (x, z, t) dz w —u x [x, 0, t) 



-ho 



where we used that u z — holds at the bottom. With the 
surface flux given by J SU rf,a; (x,t) — —Md x (S/j,), equality 
of Eq.dBl} and Eq.(|B6| implies 

d t h = d t u zlz=0 + Md 2 x {5^ z=0 ) . (B7) 

which is Eq. gDJ) at first order 0(e). 



14 



S. G. Croll, J. Appl. Polym. Sci. 23, 847 (1979). 

G. Reiter, M. Hamieh, P. Damman, S. Sclavons, [20] 
S. Gabriele, T. Vilmin, and E. Raphael, Nature Mat. 

4, 754 (2005). [21 

H. Bodiguel and C. Fretigny, Eur. Phys. J. E 19, 185 [22 
(2006). 

T. Vilmin and E. Raphael, Eur. Phys. J. E 21, 161 [23 

(2006) . 

P. Damman, S. Gabriele, S. Coppee, S. Desprez, [24 
D. Villers, T. Vilmin, E. Raphael, M. Hamieh, S. Al 
Akhrass, and G. Reiter, Phys. Rev. Lett. 99, 036101 [25 

(2007) . 

F. Ziebert and E. Raphael, Phys. Rev. E 79, 031605 [26 
(2009). 

T. Vilmin and E. Raphael, Phys. Rev. Lett. 97, 036105 [27 
(2006). 

R. Asaro and W. Tiller, Mctall. Trans. A 3, 1789 (1972). 

M. Grinfeld, Sov. Phys. Dokl. 31, 831 (1986). [28 

M. Grinfeld, J. Nonlinear Sci. 3, 35 (1993). [29 

V. Shenoy and A. Sharma, Phys. Rev. Lett. 86, 119 

(2001). 

W. Monch and S. Herminghaus, Europhys. Lett. 53, 525 [30 

(2001) . 

V. Shenoy and A. Sharma, J. Mech. Phys. Solids 50, 1155 [31 

(2002) . 

J. Sarkar, A. Sharma, and V. Shenoy, Phys. Rev. E 77, [32 
031604 (2008). [33 

D. R. Barbero and U. Steiner, Phys. Rev. Lett. 102, [34 
248303 (2009). [35 

E. Schaffer, T. Thurn-Albrecht, T. P. Russell, and 

U. Steiner, Europhys. Lett. 53, 518 (2001). [36 
P. Nozieres, in Solids far from Equilibrium, edited by [37 

C. Godreche (Cambridge University Press, Cambridge, 
1992). [38 

I. Cantat, K. Kassner, C. Misbah, and H. Miiller- [39 
Krumbhaar, Phys. Rev. E 58, 6027 (1998). 

D. Maugis, Contact, Adhesion and Rupture of Elastic 



Solids (Springer, Berlin Heidelberg, 2000). 

L. D. Landau and E. M. Lifshitz, Theory of Elasticity 

(Pergamon Press, New York, 1986). 

P. Nozieres, J. Phys. I France 3, 681 (1993). 

P. Kohlert, K. Kassner, and C. Misbah, Eur. Phys. J. B 

35, 493 (2003). 

K. Kassner, C. Misbah, J. Miiller, J. Kappey, and 

P. Kohlert, Phys. Rev. E 63, 036117 (2001). 

R. Torii and S. Balibar, J. Low Temp. Phys. 89, 391 

(1992). 

J. Berrehar, C. Caroli, C. Lapersonne-Meyer, and 

M. Schott, Phys. Rev. B 46, 13487 (1992). 

K. Sekimoto and K. Kawasaki, Physica A 154, 384 

(1989). 

C. W. Macosko, Rheology : principles, measurements, 
and applications (John Wiley & Sons Inc, New York, 
1994). 

L. H. He and L. Qiao, EPL 80, 14003 (2007). 

T. Tanaka, S. T. Sun, Y. Hirokawa, S. Katayama, 

J. Kucera, Y. Hirose, and T. Amiya, Nature 325, 796 

(1987). 

A. Oron, S. H. Davis, and S. G. Bankoff, Rev. Mod. Phys. 
69, 931 (1997). 

B. J. Spencer, P. W. Voorhees, and S. H. Davis, Phys. 
Rev. Lett. 67, 3696 (1991). 

W. M. Mullins, J. Appl. Phys. 30, 77 (1959). 

P. Miiller and A. Saul, Surf. Sci. Rep. 54, 157 (2004). 

F. Yang, J. Phys. D: Appl. Phys. 38, 3938 (2005). 

W. Song and F. Yang, J. Phys. D: Appl. Phys. 39, 4634 

(2006). 

A. Onuki, Physica A 217, 38 (1995). 

X. H. Pan, S. Q. Huang, S. W. Yu, and X. Q. Feng, J. 

Phys. D: Appl. Phys. 42, 055302 (2009). 

P. G. de Gennes, Eur. Phys. J. E 7, 31 (2002). 

Note that the case A = 1 is singluar as it yields only one 

wavenumber and additional solutions like 2sinh(.z). 



